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ABSTRACT 


In  this  stu^,  the  offset  strip  fin  array  in  a  liquid  flow  through  module  was 
examined.  The  study  used  a  model  of  the  heat  exchanger  which  is  ten  times  larger  than 
the  actual  model  used  in  electronics  cooling.  Water  and  Brayco  Micronic  889 
Polyalphaolefin  were  used  as  the  working  fluids  in  order  to  determine  the  Reynolds 
number  and  Prandtl  number  effects  on  the  Colburn  j  factor.  Also,  flow  visualization 
experiments  were  performed  to  determine  the  flow  characteristics  of  the  offset  strip  fins 
for  laminar  flow. 
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L  INTRODUCTION 


A.  ELECTRONIC  COOLING 

The  trend  in  electronic  packaging  has  been  towards  building  more  electronic 
circuits  into  one  very  small  chip.  Since  the  invention  of  the  silicon  integrated  circuit, 
there  has  been  an  increase  in  the  number  of  components  per  chip.  This  increase  in  circuit 
integration  continues  to  be  driven  by  the  need  for  increased  processing  speed.  With  each 
increase  in  the  scale  of  integration  of  microelectronic  circuits  more  effective  cooling 
systems  are  required. 

Currently,  chip  heat  fluxes  are  approaching  50W/cm^  and  heat  fluxes 
of  200  W/cm^  are  projected  [Ref  1].  Also,  it  is  known  that  chip  temperatures  must  be 
maintained  within  prescribed  operating  limits  below  85°  C.  Due  to  the  large  heat  fluxes 

and  small  areas  associated  with  current  and  projected  very  large  scale  integrated  circuits 
and  the  maxiimini  allowable  junction  temf^rature  requirement,  it  could  be  necessary  to 
cool  the  system  by  liquid  cooling  and  to  enhance  heat  transfer  by  using  extended 
surfaces. 

B.  UQUID  COOLING  METHOD 

At  lower  air  velocities,  typical  heat  flux  limits  for  packages  cooled  by 
conventional  forced  convection  are  less  than  1  W/cm^  [Ref  2].  Because  of  the  hi^  heat 
fluxes  of  today's  technology,  it  is  evident  that  air  cooling  could  not  provide  an  optimal 
solution  to  the  cooling  problem.  As  a  result  of  the  higher  density  and  thermal 
conductivity  of  liquids  in  comparison  to  air,  it  is  possible  to  enhance  heat  transfer 
coefficients  by  using  liquid  cooling. 

1.  Direct  Liquid  Cooling 

Direct  liquid  cooling  maintains  physical  contact  between  the  coolant  and  the 
electronic  components.  The  coolant  should  have  dielectric  characteristics  and  good 
chemical  compatibility  so  as  not  to  effect  the  circuit  performance. 
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exchangers.  The  cross  section  of  the  fin  is  rectangular.  The  fin  is  cut  into  small  strips  of 
length  /  in  the  flow  direction.  The  other  variables  are  the  fin  spacing  s,  fin  height  h,  and 


The  offset  strip  fins  have  high  heat  transfer  coefficients.  This  is  mainly  caused  by 
the  developing  laminar  boundary  layers.  The  interrupted  geometry  of  the  fins  prevents 
the  velocity  and  temperature  boundary  layers  from  getting  excessively  thick.  Since  the 
form  drag  increases  due  to  the  finite  thickness  of  the  fins  the  friction  factors  are  also 
high.  In  this  study  the  offset  strip  fin  array  in  a  liquid  flow  throu^  module  (LFM)  used  in 
electronics  cooling  is  examined.  The  study  uses  a  model  of  the  heat  exchanger  which  is 
ten  times  larger  than  the  actual  model.  A  typical  liquid  flow  through  module  is  shown  in 
Figure  3. 


h=  fin  height 
I  =  fin  length 
s=  fin  spacing 


Hgore  2.  Typical  Ofifeet  Strp  Fin  Geometey . 


2.  Analytical  and  Numerical  Correlations  for  Offset  Strip  Fins 

Weiting  [Ref  3]  correlated  available  experimental  heat  transfer  and  flow  friction 
data  for  22  offset  strip  fin  surfaces  using  air  as  the  working  fluid.  Weiting  defines  the 
Reynolds  number  based  on  the  hydraulic  diameter  as  follows; 
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Equation  1.1 


Here  U  is  the  mean  velocity  of  the  air  in  the  channel,  v  is  the  kinematic  viscosity 
of  the  air  and  Dh  is  the  hydraulic  diameter  {Dh  =  ^AclP)  where  Ac  represents  the 
cross-sectional  area  and  P  represents  the  perimeter  of  the  channel.  Weiting’s  correlations 

were  as  follows; 


YorReD  ^  1000, 


/■=  7.661(//Da)-® 

YotRcd  2: 2000, 

j  =  0.242(//D/,)-®  ^22(5/^^)0.089^g^.368 


Equation  1.2 
Equation  1.3 

Equation  1.4 
Equation  1.5 


Here  /  is  the  strip  length,  S  is  the  fin  thickness,  Dh  is  the  hydraulic  diameter  of  the 
passages  and  a*  is  the  ratio  of  the  width  to  height  of  the  passage.  For  Reynolds  numbers 

in  the  laminar  range  he  stated  that  the  Colburn  j  factor  was  only  a  function  of  Reynolds 
number  and  HDh  or  a*  .  He  considered  the  correlation  to  be  accurate  within  ±  10  %  , 
though  some  points  deviated  by  forty  percent. 

Manglik  and  Bergles  [Ref.  4],  gathered  all  the  previously  taken  data  for  the 
offset  strip  fins  and  developed  a  correlation  for  the  Colburn  j  factor  which  can  be  used 
for  both  the  laminar  and  turbulent  regions.  Their  definition  of  Reynolds  number  and 
hydraulic  diameter  was  the  same  widi  Weiting's  definition. 

/  =  0.6522(5//i)-«J54i(^//)0.i499(,/^)-o.0678;jg^54O3  Equation  1.6 
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In  their  correlation,  s  is  the  fin  spacing,  h  is  the  fin  height,  t  is  the  fin  thickness 

and  1  is  the  fin  length.  They  state  that  this  correlation  is  applicable  for  Prandtl  numbers  in 
the  range  of  0.5  to  15.  They  estimated  their  correlation  to  be  accurate  within  ±  20  %. 

The  correlations  for  the  transition  from  laminar  to  turbulent  region  were 

developed  by  H.  M.  Joshi  and  R.  Webb  [Ref  5].  They  used  a  numerical  solution  for  the 
Nusselt  number  and  the  friction  factor  for  a  laminar  flow  and  a  semi-emprieal  solution 


For  Ren  ^  1000, 

/  =  0.53(//Z)A)"^-^^(5//?)"®-®^i?e^'^  Equation  1.7 

For  Rec  >  Reu+  1000, 


it  is  idealized  that  the  velocity  and  temperature  boundary  layers  are  destroyed  in  the  wake 
region.  Although  this  model  is  at  best  a  first  order  approximation,  it  does  not  provide  the 
detailed  insight  into  the  systems  .  Patankar  and  Prakash  [Ref.  7],  predicted  the  j  and  f 
factors  for  the  offset  strip  fins  numerically  by  a  finite  difference  method.  They 
considered  the  laminar  boundary  layers  on  each  strip  fin.  They  first  assumed  that  strip 

fins  have  infinite  height  and  zero  thickness.  Then  they  extended  the  analysis  for  finite  fin 

thickness.  They  observed  that  at  low  Reynolds  Numbers  or  low  5/s ,  the  recirculation 
zone  behind  the  trailing  edge  was  small.  At  high  Reynolds  numbers  or  high  d/s,  this 

recirculation  zone  extended  from  the  trailing  edge  to  the  next  strip.  (See  Figure  4). 

Although  a  large  amount  of  experimentation  and  numerical  predictions  have 
been  made,  further  studies  and  improvements  are  needed  to  accurately  understand  the 
performance  of  the  offset-strip  fins. 


D,  OBJECTIVES  OF  PKESENT  STUBY 

The  first  objective  of  this  study  was  to  investigate  the  heat  transfer 
characteristics  of  a  liquid  flow  throu^  module  which  is  actually  a  ten  times  scaled  up 
model  of  the  fluid  circulation  passages  of  an  SEM-E  (Liquid  Flow  Through  Module).  By 
using  water  and  Brayco  Micronic  889  (Castrol)  Polyalphaolefin  (PAO)  as  the  working 
fluid,  experiments  were  performed  to  determine  the  Reynolds  number  and  Prandtl 
number  effects  on  the  dimensionless  heat  transfer  coefficients. 

The  second  objective  was  to  visualize  the  flow  over  the  offset  strip  fins  which 
were  used  to  enhance  the  heat  transfer.  Dye  injection  was  applied  to  observe  the  variation 
in  the  velocity  boundary  layer  across  the  fins  and  the  recirculation  of  the  flow  behind  the 
trailing  edge  of  the  fins  for  varying  flow  rate  of  the  coolant. 


8 


n.  EXPERIMENTAL  APPARATUS 


The  experimental  apparatus  included  a  system  to  circulate  the  coolant  through 
the  test  section  and  another  system  to  collect  the  necessary  data  to  be  able  to  calculate  the 
heat  transfer  and  flow  characteristics  of  the  test  section. 

A.  COOLANT  CIRCULATING  SYSTEM 

A  schematic  of  the  coolant  circulation  system  is  shown  in  Figure  5. 


Coolant  Circulating  System . 


1.  Liquid  Tank 

In  this  study  first  water  and  then  Brayco  Micronic  889  (Castrol)  Polyalphaolefin 
(PAO)  were  used  as  the  working  fluids.  In  both  case,  the  coolant  was  provided  by  an 


The  flowmeter  was  calibrated  for  both  of  the  coolants.  The  calibration  of  the 
flowmeter  for  water  was  made  in  two  different  ways.  The  first  one  consisted  of  a 
comparison  of  the  average  voltage  output  to  the  average  time  to  collect  a  specific  volume 
of  water.  The  second  method  was  based  on  the  collected  mass  of  the  water  in  a  specific 
time  for  varying  pump  speeds.  Even  though  the  second  method  produced  a  more  linear 
calibration  curve  than  the  first  one,  the  flow  rate  measurements,  especially  in  the  range  of 
middle  Reynolds  numbers  (200-600),  were  very  close  to  each  other.  The  calibration  of 
the  flowmeter  for  the  PAO  was  based  on  the  first  method.  All  flowmeter  calibration 
curves  were  presented  as  Appendix  A. 

4.  Heat  Exchanger 

The  flowmeter  outlet  was  connected  to  a  manifold  and  divided  into  three  flow 
paths  before  entering  the  offset  strip  fin  heat  exchanger. 

a.  Offsa  IStr^  Fins 

The  test  section  which  includes  the  base  plate  and  the  oflfeet  strip  fin 
assembly  was  first  constructed  by  Lt.  Jeffrey  Masterson  [Ref  8].  A  picture  of  the  test 
section  is  shown  in  Figure  6.  The  fin  assembly  was  a  10  times  scaled  up  model  of  the 
horizontal  passages  of  a  liquid  flow  through  module  used  in  electronics  cooling.  The  base 
plate  consisted  of  three  parts.  (The  inlet  plenum,  outlet  plenum  and  the  middle  part  with 
fins).  The  inlet  and  outlet  plenum  were  made  fi-om  two  pieces  of  a  12.7  mm  aluminum 
6061  alloy  plate.  They  were  milled  down  to  the  base  plate  thickness  of  10.16  mm  The 
middle  part  of  the  base  plate,  which  included  the  offset  strip  fins,  was  also  made  fi-om 
aluminum  6061  alloy.  Instead  of  welding  the  fins  onto  the  base  plate,  the  base  and  the 
fins  were  milled  fi^om  a  solid  piece,  not  to  leave  weld  deposits  and  affect  the  fin 
geometry.  A  2.54  cm  thick  aluminum  plate  was  cut  into  31.75  mm  wide  strips.  These 
strips  were  then  milled  to  the  width  of  the  plate.  The  fins  were  cut  15.24  mm  deep  into 
the  aluminum  and  15  rows  were  obtained.  Each  row  had  19  fins  where  the  spacing 
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Figure  6.  Test  Section. 
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between  the  fins  was  11.76  mm.  A  6.64  mm  offset  was  allowed  on  one  end  that  when 
alternated  created  the  offset  pattern  of  succeeding  fin  rows.  On  each  end  of  the  base  plate 
an  additional  12.7  mm  offset  of  base  plate  thickness  was  milled  to  attach  the  plexiglas 
cover  [Ref.  8].  A  schematic  of  the  finned  section  is  shown  in  Figure  7. 


figure  7.  Finned  Section  Geometry. 


The  inlet  plenum,  firmed  section  and  the  outlet  plenum  were  connected  to  each 
other  by  means  of  two  threaded  rods.  The  holes  to  accommodate  the  threaded  rods  were 
drilled  in  the  lengthwise  direction  through  each  fin  row  and  plenum  piece.  Rubber  gasket 
was  placed  between  the  first  row  and  the  inlet  plenum  sections  and  between  the  last  row 
and  the  outlet  plenum  sections  to  reduce  the  longitudinal  conduction  losses  fi’om  the 
finned  section  to  the  plenum.  The  final  assembly  was  held  and  tightened  by  using  nuts  at 
each  end  of  the  rods  [Ref  8]. 


axial  ends  of  the  plexiglas  cover.  Inlets  and  outlets  consisted  of  3.175  mm  pipe  threads 
drilled  and  tapped  into  the  plexiglas  wall  with  6.35  mm  tygon  adaptors  on  the  opposite 
end  of  the  wall.  One  inlet  and  one  outlet  were  aligned  on  the  longitudinal  axis  of  the  test 
plate.  The  two  other  inlets  and  outlets  were  placed  37.9  mm  from  the  outer  edges  of  the 
plate.  Valves  were  placed  on  the  two  other  outlets  so  they  could  be  closed  when 
operating  the  test  section  as  needed  or  opened  to  drain  the  system  after  each  experiment. 
During  the  first  runs  these  outlets  with  valves  were  kept  closed  and  the  inlet  to  the 
circulating  bath  was  provided  by  only  one  outlet  of  the  heat  exchanger.  But  since  the 
pressure  in  the  heat  exchanger  was  high  especially  for  Re  >  500,  one  of  the  outlet  valves 
was  also  kept  open  vsdiich  provided  the  inlet  to  the  circulating  bath  from  two  flow  paths. 
This  decreased  the  hig^  back  pressure  in  the  heat  exchanger. 

Another  step  before  attaching  the  cover  to  the  base  plate  was  to  drill  the 
dye  injection  and  temperature  probe  points  into  the  plexiglas  cover.  The  details  about  the 
dye  injection  points  will  be  given  in  the  next  section. 

Four  temperature  probe  holes  were  placed  above  both  the  inlet  and  outlet 
plenum  regions  to  measure  the  inlet  and  outlet  temperatures  of  the  coolant.  These  holes 
were  very  tiny  and  filled  with  silicone  before  feeding  the  thermocouples  through  the 
holes. 

Prior  to  attaching  the  cover  to  the  base  a  flow  straightener  was  also  fixed 
between  the  inlet  temperature  probe  and  dye  injection  holes.  The  original  flow 
straightener  which  was  devised  using  38. 1  mm  long  straws  and  a  plastic  mesh  was  glued 
to  the  plexiglas  cover.  The  purpose  was  to  provide  a  uniform  flow  velocity  profile  before 
the  first  fin  row.  The  distance  between  the  first  fin  row  and  the  edge  of  the  flow 
straightener  was  about  28  mm.  After  drilling  all  the  necessary  holes  and  fixing  the  flow 
straightener,  the  sides  were  attached  to  the  top  piece  using  acrylic  cement  and  secured 
with  screws.  During  the  first  runs  leakage  was  noted  from  the  temperature  probe  holes 
and  the  side  walls  of  the  plexiglas  cover.  Silicone  was  run  along  the  sides  and  refilled 
into  the  temperature  probe  holes.  When  the  heat  transfer  experiments  were  completed  a 
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5ig5sr@9.  Dye  Lyection  System. 
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Five  dye  injection  points  were  placed  in  the  plexiglas  cover.  They  were  1.6  mm 
holes  that  penetrated  the  plexiglas  cover  just  before  the  first  row  fins  and  after  flow 
straightener  assembly.  Five  adaptors  were  placed  into  these  holes.  Inside  the  plexiglas 
cover,  stainless  steel  surgical  tubing  was  passed  through  the  adaptors.  This  tubing  was 
bent  into  a  90  degree  bend  and  the  ends  were  smoothed  and  shaped  to  a  cone  in  order  to 
direct  the  dye  flow  between  the  fins.  The  total  length  of  a  steel  tubing  was  about  50  mm. 
The  other  ends  of  the  adaptors  which  were  left  on  the  outside  of  the  plexiglas  cover  were 
connected  to  a  manifold  by  using  190  mm  long  plastic  tubes.  The  manifold  was  used  to 
distribute  the  dye  in  the  reservoir  to  the  desired  injection  tubings.  The  reservoir  was  a 
small  plastic  bottle  and  was  held  about  1  meter  above  the  test  section  to  provide  a 
positive  pressure.  A  small  syringe  also  served  as  an  ink  reservoir. 

C.  DATA  COLLECTING  SYSTEM 

A  heater  was  used  to  heat  the  test  section  and  the  temperatures  were  measured 
by  means  of  thermocouples.  Thermocouples  were  connected  to  a  data  acquisition  system 
which  also  measured  the  required  voltages.  Flow  visualization  data  was  recorded  by 
using  a  35  mm  camera. 

1.  Heater 

Printed  circuit  boards  usually  dissipate  constant  heat  flux.  This  was  simulated 
by  using  an  electrical  heater  which  was  a  Minco  foil  backed  patch.  The  heater  pad  was 
25.4  cm  by  30.48  cm,  resulting  in  a  heater  surface  area  of 774.2  cm^ . 

The  heater  pad  was  fixed  to  the  bottom  of  the  aluminum  plate  with  pressure 
sensitive  adhesive.  The  pad  covered  the  finned  area  section  except  3.24  mm  on  the  sides 

and  9.525  mm  at  each  end.  A  Kepco  0-100  V,  0-5  A  power  supply  provided  DC  power  to 
the  heater.  A  250  Watt,  2±0.01  n  precision  resistor  was  connected  in  series  with  the 
heater.  The  purpose  of  the  precision  resistor  was  to  determine  the  power  to  the  heater. 
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Equation  2. 1 


h  =  Ip  Equation  2.2 

Powers hxVh  Equation  2.3 

2.  Thermocooples 

The  temperature  of  the  aluminum  base  plate  was  measured  by  using  23  Copper 
Constantan  (T-Type)  cement-on  thermocouples.  They  were  placed  on  the  bottom  siuface 
of  the  finned  base  plate  as  seen  in  Figure  11.  Since  the  thermocouples  were  placed 
between  the  heater  pad  and  the  base  plate,  flat  ribbon  thermocouples  were  chosen  to 
provide  better  surface  contact.  The  thermocouples  were  glued  using  Omegabond  101 
which  was  a  high  thermal  conductivity  epoxy. 

Eight  Copper  Constantan  (T-Type)  thermocouples  were  fed  through  the  holes 
drilled  into  the  plexiglas  cover  to  measure  the  inlet  and  outlet  coolant  temperatures.  Two 
thermocouples  were  affixed  between  the  plexiglas  cover  and  the  insulation  layer  where 
the  two  others  were  moimted  on  top  of  the  insulation  layer  to  be  able  to  estimate  the  heat 
losses  through  the  cover  and  the  insulation. 

Another  Copper  Constantan  (T-Type)  thermocouple  was  placed  into  the 
refrigerated  bath  to  measure  the  coolant  temperature  in  the  bath.  All  the  thermocouples 
were  then  connected  to  the  Data  Acquisition  System. 

3.  Data  Acquisition  System 

All  the  thermocouples  and  voltages  were  measured  by  an  Hewlett  Packard  3852 
Data  Acquisition/Control  Unit.  Two  24-channel,  thermocouple  compensated  high  speed 
multiplexed  boards  (HP  44713  A)  measured  a  total  number  of  34  thermocouple  junction 
temperatures  and  a  24-channel  high  speed  multiplexed  board  (HP  44705  A)  measured  the 
voltage  drops  across  the  flowmeter,  precision  resistor  and  the  power  supply. 


All  thermocouple  extensions  were  labeled  and  plugged  into  channel  numbers  0 
through  17  of  both  HP  44713  A  boards  and  the  voltage  drops  were  measured  by  the 
channels  0  through  2  of  the  HP  44705  A  board. 

The  Data  Acquisition  Unit  was  controlled  by  a  Personal  Computer.  An  interface 
board,  G  PIB-PCIIA  IEEE-488,  manufactured  by  National  Instruments  was  installed  into 
this  computer.  A  program  attached  as  Appendix  B  was  written  in  QBASIC  code  to  store 
the  data  for  later  computation  and  display.  A  schematic  of  the  entire  experimental 
apparatus  is  presented  in  Figure  12. 
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m.  EXPERIMENTAL  PROCEDURE 


A.  HEAT  TRANSFER  EXPERIMENTS 

Heat  transfer  experiments  were  performed  for  both  water  and  Brayco  Micronic 
889  (Castrol)  Polyalphaolefin  (PAO).  Since  the  purpose  of  these  experiments  was  to 
determine  the  Colburn  j  factor  for  varying  Reynolds  numbers,  the  flowrate  was  varied  to 
be  able  to  obtain  the  desired  Reynolds  numbers  of  interest. 

The  flow  velocity  (Reynolds  number)  was  limited  by  the  maximnin  pressure  the 
heat  exchanger  could  withstand.  With  the  power  to  the  heater  held  constant,  the  pump 
speed  was  varied  between  10-70  %  of  its  maximum  value  which  corresponded  to 
Reynolds  numbers  in  the  range  of  50-650  for  water  and  Reynolds  numbers  in  the  range  of 
10-80  for  Brayco  Micronic  889  Polyalphaolefin.  This  procedure  was  repeated  for  various 
power  inputs  to  the  heater  (100-250  Watts).  The  power  input  to  the  heater  was  limited  to 
250  Watts  due  to  precision  resistor  limitations. 

The  bath  temperature  was  maintained  at  around  20"C.  As  stated  previously, 
^ftcr  the  system  reached  the  steady  state  condition,  the  necessary  data  was  collected  by 

the  Data  Acquisition  Unit  and  used  to  determine  the  heat  transfer  characteristics  of  the 
test  section.  For  the  calculations  the  following  steps  were  followed: 

1.  Reynolds  Number  Calculations 

All  the  calculations  were  based  on  the  hydraulic  diameter.  The  hydraulic 
diameter  relates  the  cross-sectional  area  of  a  channel  to  its  heat  transfer  area.  For  the 
present  study  different  area  defimtions  were  defined  to  determine  the  hydraulic  diameter 
and  the  Reynolds  number. 

•  Channel  Cross  Sectional  Area  (Ac): 

The  product  of  the  fin  height  h,  and  the  fin  spacing  s,  gives  the  cross-sectional 
area  of  each  channel. 
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Equation  3.4 


AH  represents  the  characteristics  length  of  each  channel.  As  stated  previously, 
the  volumetric  flow  rate,  Q  was  determined  by  using  the  flowmeter  calibration  curves 

obtained  for  both  water  and  Brayco  Micronic  889  Polyalphaolefin.  (See  Appendix  A). 
Dividing  the  volumetric  flow  rate  by  the  total  cross-sectional  area  A  /,  the  velocity  V 
through  the  channels  was  determined. 


Equation  3.5 


The  fluid  properties  were  determined  by  using  the  film  temperature,  Tf.  The 
density  p,  viscosity  |A,  specific  heat  Cp,  thermal  conductivity  Jt/  and  Prandtl  number  Pr 
of  water  taken  from  Ref  9  were  curve  fitted  in  the  range  of  0-80  °C.The  same  procedure 

was  repeated  for  Brayco  Micronic  889  Polyalphaolefin.  The  properties  of  this  coolant  are 
tabulated  as  Appendix  C  according  to  the  data  provided  by  the  manufacturer,  Castrol 
[Ref.  10]. 

Finally  the  Reynolds  number  based  on  the  hydraulic  diameter  was  determined 
as  follows; 


Rei)  = 


pxDfjxV 

ii 


Equation  3.6 


2.  Colburn  j  Factor  Calcnlations 

The  film  temperature  Ti,  was  averaged  between  the  plate  average  temperature 
Tavg,  and  the  coolant  inlet  temperature  Tm.  As  explained  in  the  section  about 

thermocouples,  the  inlet  temperature  of  the  coolant  was  measured  at  four  locations.  The 

plate  average  temperature  was  the  average  temperature  of  the  23  thermocouples  which 

were  placed  between  the  heater  pad  and  the  base  plate. 


T avg-^Tff, 
2 


Equation  3.7 


rr,  T412+T413+7’414+7415  _ 

im-  4  Equation  3.8 
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7’(100-111)+7’(40(>-410) 
Tavg=  23 


Equation  3.9 


(The  subscript  numbers  of  T  in  Equations  3.8  and  3.9  represent  the  thermocouple 
numbers). 

The  electrical  power  input  PeUc  to  the  heater  was  calculated  by  using  the 
equations  2.1, 2.2  and  2.3.  The  temperature  difference  through  the  insulation  layer  AT;^, 

used  to  estimate  the  heat  losses  from  the  heat  exchanger. 

T4i6+r4i7-7’n6-rii7  ^ 

^ - —  Equation  3.10 


Since  the  thickness  of  the  insulation  tmsui  was  about  0.01905  m  and  the  thermal 
conductivity  kj^ui  was  0.04  W/mK.  The  lost  or  gained  heat  was  determined  by  using  the 

following  formula. 


AP  = 


kjnsul. 


Equation  3.11 


Here  Ah  represents  the  area  covered  with  the  insulation  layer.  The  actual  power 
F,  to  the  system  was  obtained  by  adding  AF  to  the  electrical  power. 


F  =  Fefec,+AF 


Equation  3.12 


The  next  step  was  to  calculate  the  heat  flux.  Since  the  heater  dimensions  and 

the  power  input  to  the  heater  was  held  constant  for  a  set  of  experiments,  the  heat  flux  q” 
was  also  constant.  The  heat  flux  and  the  average  Nusselt  number  Nuavg  were  calculated 

as  follows: 

p 

^”~Ah  Equation  3.13 
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q”^Dh 


Equation  3.14 


Finally,  the  experimental  Colburn  j  factor  was  calculated  to  relate  the 
dimensionless  heat  transfer  coefficient  to  the  Reynolds  number. 

NUavg  _ 

B.  FLOW  VISUALIZATION  EXPERIMENTS 

The  procedure  for  these  experiments  was  to  change  the  flow  velocity  of  the 
water  and  visualize  the  flow  on  the  offset  strip  fins.  Only  water  was  used  as  the  working 
fluid  and  the  power  input  to  the  heater  was  set  to  zero  during  the  visualization 
experiments.  The  pump  speed  was  varied  to  obtain  the  same  Reynolds  numbers  of  the 
heat  transfer  experiments  performed  with  water  (Re  50-650).  Each  fin  row  was  labeled 
on  the  plexiglass  cover  to  be  able  to  determine  the  location  of  the  changes  in  the  velocity 
profile.  First  blue  ink  and  then  white  dye  were  injected  into  the  plexiglas  cover  and 
directed  through  the  fin  array.  Photographs  were  taken  for  each  Reynolds  number 
settings.  Several  photographs  should  be  taken  to  be  able  to  observe  all  the  fin  rows. 


IV.  RESULTS 


A.  HEAT  TRANSFER  RESULTS 

When  water  was  used  as  the  working  coolant  the  power  supply  was  set  to  100, 
150,  200  and  250  Watts.  The  data  obtained  and  the  calculated  Colburn  j  factors  are 
tabulated  for  each  power  setting  in  Tables  14.  Also  some  of  the  experimental  data  were 
plotted  for  varying  Reynolds  numbers  and  power  inputs. 

The  differences  between  the  water  inlet,  and  water  outlet  temperatures  are 
plotted  on  a  graph  (Figure  13),  while  the  differences  between  the  average  plate  and  water 
inlet  temperatures  are  plotted  on  another  graph  (Figure  14),  for  all  power  settings.  Figure 
15,  shows  the  plate  average  temperatures  for  varying  Reynolds  numbers.  The  calculated 
Colburn  j  factors  are  also  plotted  separately  for  each  power  setting  in  Figures  16-19. 

The  same  procedure  was  repeated  for  the  coolant  Brayco  Micronic  889 
Polyalphaolefin.  This  time  the  power  supply  was  set  to  100, 175  and  250  Watts.  The  data 
obtained  were  tabulated  in  Tables  5-7.  Temperature  differences  mentioned  above  are 
plotted  in  Figures  20  and  21.  The  calculated  Colburn  j  factors  for  the  coolant  PAO  are 
shown  in  Figures  22-24.  The  correlations  developed  by  Weiting  [Ref  3],  and  Manglik 
and  Bergles  [Ref  4],  were  used  to  predict  the  Colburn  j  factor  for  the  present  study.  (See 
Chapter  I,  Section  C).  These  predicted  Colburn  j  factors  were  plotted  in  Figures  25-28  to 
be  able  to  compare  with  the  experimental  results  obtained  for  various  power  inputs. 
These  figures  also  included  the  results  for  the  coolant  PAO  to  determine  the  Prandtl 
number  effect  on  the  Colburn  j  factor.  It  was  noted  that  the  Prandtl  number  of  the  PAO 
(Brayco  Micronic  889)  was  around  230-250  and  the  Prandtl  number  of  water  was  around 
5-7  during  the  experiments. 

The  data  and  the  graphs  showed  that  as  the  Reynolds  number  increased,  the 
Nusselt  number  also  increased  for  both  of  the  coolants.  It  was  noted  that  higher  power 
inputs  produced  higher  Nusselt  numbers.  For  the  same  power  input  and  Reynolds 
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accuracy  zone.  The  difference  between  the  predicted  and  calculated  Colburn  j  factors  for 
the  power  input  of 250  Watts  were  tabulated  in  Table  8. 

Since  the  predictive  correlations  were  applicable  for  Prandtl  numbers  in  the 
range  of  0.5  to  15,  and  the  Prandtl  number  of  Brayco  Micronic  889  (PAO),  was  about 
230-250  during  the  experiments,  there  was  not  a  close  agreement  between  the  predicted 
and  experimental  results  when  this  coolant  (PAO)  was  used  as  the  working  fluid.  Even 
though  the  experimental  results  were  higher  than  the  predicted  values  the  curves  obtained 
for  both  of  the  experimental  and  predicted  Colburn  j  factors  showed  a  similar  tendency. 
(See  Figures  25-28).  This  validated  the  Prandtl  number  limitations  of  these  predicted 
correlations. 

Finally,  the  Colburn  j  factors  of  water  and  Brayco  Micronic  889  (PAO)  were 
compared  to  each  other.  When  the  power  input  was  250  Watts  and  the  Reynolds  number 
was  around  57,  the  calculated  Colburn  j  factor  was  0.0856  for  water  (Pr  =  6.03),  and 
0.3280  for  the  PAO  (Pr  =  238.72),  it  was  noted  that  the  ratio  between  the  Colburn  j 
factors  was  about  the  same  to  the  1/3  to  the  power  of  the  ratio  between  the  Prandtl 
numbers. 


jpAO  _  0.3280  __ 
jwater  0.0856 


The  difference  between  the  ratios  was  about  11%. 
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Power  Input »250  Watts  <Heat  Flux  ®  3230  W/pi® 
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.  Experimental  Data  for  Water  (Power  Input  =  250  W). 


Power  Input =175  Watts  <Heat  Flux  ®  2260  W/n* 


Fluid:Water- 
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Difference  vs.  Reynolds  Number. 


100  Uatts 
150  Uatts 


.  Colburn  j  Factor  vs.  Re3?nolds  Ni 


0.09  _ _ Colburn  j 


r  vs.  Re5?nolds  Number  (Power  =  2(X)W). 


=  Temperal 


Working  Fluid:  Brayco  Micronic  889 
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.  Colburn  j  Factor  vs.  Reynolds  Number  (Power  =103  W). 


Colburn  j  Factor  vs.  Roynoids  Number 


Reynolds  Number 
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V.  CONCLUSIONS 


The  heat  transfer  experiments  of  this  study  showed  that,  when  the  working  fluid 
was  water,  the  experimental  values  come  to  close  agreement  to  the  predicted  values.  But 
when  Brayco  Micronic  889  (PAO)  was  used  as  the  working  fluid,  the  experimental 
results  were  significantly  higher  than  the  predicted  results.  This  provides  evidence  that 
the  previous  correlations  were  valid  for  water  but  not  for  PAO. 

The  comparison  of  the  predicted  dimensionless  heat  transfer  coefficients  to  the 
experimental  values  validated  the  Prandtl  number  dependency  of  the  Colburn  j  factor.  It 
was  obtained  that  the  increase  in  the  Colburn  j  factor  was  about  the  same  with  1/3  to  the 
power  of  the  increase  in  the  Prandtl  number  of  the  working  fluid.  This  also  proved  the 
superiority  of  the  higher  Prandtl  number  fluids  on  the  lower  Prandtl  number  fluids  to 
increase  the  performance  of  an  heat  exchanger. 

Flow  visualization  experiments  showed  that,  the  interrupted  geometry  of  the 
fins  prevented  the  velocity  boundary  layers  from  getting  excessively  thick.  Observing 
these  laminar  boundary  layers  along  the  offset  strip  fins  helped  to  explain  the 
enhancement  in  the  heat  transfer  of  the  system  for  varying  flow  rates. 

In  general,  this  study  will  help  the  researchers  to  predict  the  dimensionless  heat 
transfer  coefficient  of  a  ten  times  scaled  up  model  of  a  liquid  flow  through  module  used 
in  electronics  cooling. 
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VL  RECOMMENDATIONS 


In  any  continuation  of  this  study  the  followings  are  recommended; 

•  Clean  the  water  deposits  on  the  base  plate  and  offset  strip  fins  and  prepare  the 
test  section  for  the  new  study. 

•  Redesign  and  manufacture  a  stronger,  pressure  resistant  plexiglas  cover  for  the 
heat  exchanger. 

•  Insert  additional  thermocouples  into  the  inlet  and  outlet  of  the  test  section. 
Remove  the  heater  and  calibrate  the  thermocouples  mounted  on  the  bottom  of  the  base 

plate  fin  assembly  to  a  zero  reference  point. 

•  Repeat  the  experiments  of  the  present  study  by  using  a  pump  with  higher 
capacity  to  investigate  the  heat  transfer  and  flow  characteristics  of  the  system  for  the 

turbulent  flow. 

•  Develop  a  numerical  solution  to  the  flow  through  the  offset  strip  fins  and 
compare  the  solution  with  the  experimental  results. 

•  Investigate  the  Fanning  friction  factor  and  the  effect  of  the  fin  thickness  on  the 
form  drag. 
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APPENDIX  A.  FLOWMETER  CALIBRATION 


The  flowmeter  was  calibrated  in  two  different  ways.  The  first  one  consisted  of  a 
comparison  of  the  average  voltage  output  to  the  average  time  to  collect  a  specific 
volume  of  water  as  shown  in  Figure  3 1 . 


Th©  second  msthod  was  bas©d  on  the  collected  mass  of  the  water  in  a  specific 
time  for  varying  pump  speeds  as  shown  in  Figure  32. 


Volumetric  Flowrate  (Umln) 


The  calibration  of  the  flowmeter  for  Brayco  Micronic  889  Polyalphaolefin  (PAO) 
was  based  on  the  first  method  which  consisted  the  comparison  of  the  flow  percentage 
readings  to  the  average  time  to  collect  a  specific  volume  of  PAO  as  shown  in  Figure  33. 
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APPENDIX  B.  QBASIC  PROGRAM  CODE 


CLS 

COLOR  14, 9 

PBONT  "  Program  for  National  Instruments  PC2  IEEE-488  BOARD  TO  HP-3852" 
PRINT "  This  uses  the  Universal  Language  Interface" 

PRINT "  ULLCOM  must  be  run  prior  to  running  the  program" 

PRINT "  The  program  name  is  COLBURNJ.BAS" 

PRINT  "  Modified  by  Ata  KOMRAL" 

OPEN  "GPIBO"  FOR  OUTPUT  AS  #1 
OPEN  "GPIBO"  FOR  INPUT  AS  #2 
PRINT  #1,  "ABORT" 

PRINT  #1,  "RESET" 

PRINT  #1,  "GPIBEOS  CR  LF" 

PRINT  #1,  "CLEAR  9"  STARTIT; 

DIM  TEMP(20),  TMP(20),  V(3),  TIN(33),  TOUT(33),  DIFF(33),  TAV(33),  FILM(33), 
TINSULA(33),  TPGLSS(33) 

PRINT  #1,  "REMOTE  9" 

PRINT  #1,  "OUTPUT  9;RESET  700" 

PRINT  #1,  "OUTPUT  9;USE  700" 

INPUT  "ENTER  PUMP  SETTING ",  SET 
TOTAV  =  0! 

TOTIN  =  0! 

TOTOUT  =  0! 

TOTDIFF  =  0! 

TOTFILM  =  0! 

TOTCVR  =  0! 

TOTINSL  =  0! 

N  =  33 

PRINT  TIMES; 

FORJ=lTON 

PRINT  #1,  "OUTPUT  9;CONF  TEMPT  " 

FORI  =  0TO17 
CHNL$  =  STR$(I+  100) 

MSGS  =  "OUTPUT  9;MEAS  TEMPT"  +  CHNLS 
PRINT  #1,  MSGS 
PRINT  #1,  "ENTER  9" 

INPUT  #2,  IMPS 
TMP(I)  =  VAL(TMPS) 

NEXTI 

PRINT  #1,  "OUTPUT  9;  CONF  TEMPT" 

FORI  =  0TO17 


CHN$  =  STR$(I  +  400) 

MSGE$  =  "OUTPUT  9;  MEAS  TEMPT"  +  CHN$ 
PRINT  #1,MSGE$ 

PRINT  #1,  "ENTER  9" 

INPUT  #2,  TEMPS 
TEMP(I)  =  VAL(TEMP$) 

NEXTI 

PRINT  #1,  "OUTPUT  9;CONF  DCV" 

PRINT  #1,  "OUTPUT  9;USE  700" 

FORI  =  0TO2 
CH$  =  STR$(I+500) 

MESSS  =  "OUTPUT  9;MEAS  DCV"  +  CHS 
PRINT  #1,MESS$ 

PRINT  #1,  "ENTER  9" 

INPUT  #2,  VOLTS 
V(I)  =  VAL(VOLT$) 

NEXTI 

PRINT#!,  "CLEAR  9" 

PRINT  #1,  "LOCAL  9" 

TPLATEI  =  0! 

1  =  0 

DOWfflLEI<12 
1  =  1+1 

TPLATEI  =  TPLATEI  +  TMP(I  - 1) 

LOOP 

TPLATE2  =  0! 

1  =  0 

DOWHILEI<ll 
1  =  1+  1 

TPLATE2  =  TPLATE2  +  TEMP(I  - 1) 

LOOP 

TSURF  =  (TPLATEI  +  TPLATE2)  /  (12  +  11) 

TIN  =  0! 

1  =  0 

DOWHILEI<4 
1  =  1+  1 

TIN  =  TIN  +  TEMP(I+  11) 

LOOP 

INTMP  =  TIN/4 
TOUT  =  0! 

1  =  0 

DOWHILEI<4 
1  =  1+  1 
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TOUT  =  TOUT  +  TMP(I  +  1 1) 

LOOP 

OUTTMP  =  TOUT/4 
TPGLASS  =  0! 

1  =  0 

DOWHILEI<2 
1  =  1+1 

TPGLASS  =  TPGLASS  +  TMP(I  +  15) 
LOOP 

TCVR  =  TPGLASS/ 2 
TINS  =  0! 

1  =  0 

DOWHILEI<2 
1  =  1+  1 

TINS  =  TINS  +  TEMP(I  +  15) 

LOOP 

TINSUL  =  TINS/2 
TAV(J)  =  TSURF 
TIN(J)  =  INTMP 
TOUT(J)  =  OUTTMP 
DIFF(J)  =  OUTTMP  -  INTMP 
FILM(J)  =  (TSURF  +  INTMP)  /  2 
TPGLSS(J)  =  TCVR 
TINSULA(J)  =  TINSUL 
TOTAV  =  TOTAV  +  TAV(J) 

TOTIN  =  TOTIN  +  TIN(J) 

TOTOUT  =  TOTOUT  +  TOUT(J) 
TOTDIFF  =  TOTDIFF  +  DIFF(J) 

TOTFILM  =  TOTFILM  +  FILM(J) 
TOTCVR  =  TOTCVR  +  TPGLSS(J) 
TOTINSL  =  TOTINSL  +  TINSULA(J) 
NEXTJ 

TSURFACE  =  TOTAV  /  N 

INLET  =  TOTIN /N 

OUTLET  =  TOTOUT  /  N 

DELTA  =  TOTDIFF /N 

FILMTEMP  =  TOTFILM  /N 

TCOVER  =  TOTCVR  /  N 

TINSULATION  =  TOTINSL  /  N 

DELTINSUL  =  TINSULATION  -  TCOVER 

TEMDff  =  TSURFACE  -  INLET 

R  =  2 

P  =  V(2)*(V(1)/R) 
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FOR  I  =  0  TO  2 

PRINT  USING  ’•###.##";  V(I); 

NEXT  I 

PRINT  TWR=’’ 

PRINT  USING  "##.##";? 

FORI=0TO3 
PRINT  I  +  100; 

PRINT  USING  "###.###’’;  TMP(I); 

PRINT  I  +  104; 

PRINT  USING  "###.###’’;  TMP(I  +  4); 

PRINT  I  + 108; 

PRINT  USING  "###.###";  TMP(I  +  8); 

PRINT  I +112; 

PRINT  USING  "###.###";  TMP(I  +  12) 

NEXT  I 
FORI  =  0TO3 
PRINT  I +  400; 

PRINT  USING  "###.###";  TEMP(I); 

PRINT  I +  404; 

PRINT  USING  "###.###";  TEMP(I  +  4); 

PRINT  I +  408; 

PRINT  USING  "###.###";  TEMP(I  +  8); 

PRINT  I +  412; 

PRINT  USING  "###.###";  TEMP(I  +  12) 

NEXT  I 

PRINT  #1,  "OUTPUT  93EEP" 

PRINT#1,"(XEAR9" 

PRINT#l,’TOCAL9" 

z=19.5 

1  =  . 03175 

s=.01175 

t  =  . 00152 

h  =  . 01524 

Ah  =  .07742 

Ac  =  s  h 

Af  =  Ac  z 

At  =  s*I  +  2*h*l  +  2*t*h 
Dh  =  (4«Ac)/(At/l) 

vis  =  -(3.7292  / 10  ^  12)  ^  FILMTEMP  ^  3  +  (7.125  /  10  ^  10)  *  FILMTEMP  2 
-  (5.0808  /  10  ^  8)  *  FILMTEMP  +  (1.7758  / 10  ^  6) 

kf  =  (8.3333  / 10  ^  8)  *  FILMTEMP  ^  3  -  (2.1429  / 10  ^  5)  *  FILMTEMP  ^  2 
+  (2.631  / 10  ^  3)  *  FILMTEMP  +  .55206 

Pr  =  (4.4271  / 10  ^  7)  FILMTEMP  ^  4  -  (1 .0604  / 10  ^  4)  FE.MTEMP  3 
+  (9.9979  / 10  ^  3)  FILMTEMP  ^  2  -  (4.9008  / 10)  ♦  FILMTEMP  +  13.6 
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den  =  (2/10^5)*  FILMTEMP  ^  3  -  (6.407  / 10  ^  3)  *  FILMTEMP  ^  2 
+  (3.207  / 10  ^  2)  ♦  FILMTEMP  +  1002.3 

Cp  =  (5.2604  /  10  ^  6)  ♦  FILMTEMP  ^  4  -  ( 1 . 1 1 67  / 1 0  ^  3)  *  FILMTEMP ""  3 

+  (9.3021  / 10  ^  2)  *  FILMTEMP  ^  2  -  (3.2558)  *  FILMTEMP  +  4217.8 

m  =  . 001 8*  SET -.0073 

Q  =  m  /  den 

VEL  =  Q/Af 

Re  =  (VEL*Dh)/vis 

kins  =  .04 

thick  =  .75  *2.54/ 100 

DELP  =  kins  ♦  Ah  ♦  DELTINSUL  /  thick 

PTOT  =  P  +  DELP 

FLUX-P/Ah 

Nu  =  (FLUX  *  Dh)/kf/  TEMDIF 
Je  =  Nu/Re/(Pr^(l/3)) 

PRINT  TIMES 
PRINT  "TAVG="; 

PRINT  USING  "###.###";  TSURFACE; 

PRINT  "TIN="; 

PRINT  USING  ”###.###";  INLET; 

PRINT  "TOUT="; 

PRINT  USING  "###.###";  OUTLET; 

PRINT  "DELTA="; 

PRINT  USING  "###.###";  DELTA; 

PRINT  "TFILM="; 

PRINT  USING "  ###.###";  FILMTEMP; 

PRINT  "TDIF="; 

PRINT  USING "  ###.###”;  TEMDIF; 

PRINT  "TCOVER="; 

PRINT  USING  ’’  ###.###";  TCOVER; 

PRINT  "TINSUL="; 

PRINT  USING "  ###.###";  TINSULATION; 

PRINT  "visc="; 

PRINT  USING  "###.######";  vis; 

PRINT  "cond="; 

PRINT  USING  ”##.######";  kf; 

PRINT  "Prandtl="; 

PRINT  USING  "##.######";  Pr, 

PRINT  "Density="; 

PRINT  USING  "####.#####";  den; 

PRINT  ’’Cp="; 

PRINT  USING  '•####.####";  Cp; 

PRINT  TOWER=’'; 
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PRINT  USING  ''####.####••;  PTOT; 

PRINT  "m=”; 

PRINT  USING  "###.####";  m; 

PRINT  "VEL="; 

PRINT  USING  "###.#####";  VEL; 

PRINT  "Re="; 

PRINT  USING  "###.#####";  Re; 

PRINT  "FLUX="; 

PRINT  USING  "####.#####";  FLUX; 

PRINT  "Nu="; 

PRINT  USING  "####.#####";  Nu; 

PRINT 

PRINT  USING  "####.#####";  Je; 

INPUT  "  DO  YOU  WISH  TO  SAVE  DATA  ?  ",  SAVTTS 
IF  SAVrrS  =  "n"  THEN  GOTO  QUES 
IF  SAVrrS  =  •'N"  THEN  GOTO  QUES 
FILEIT: 

PRINT "  IF  DATA  FILE  EXISTS  OLD  DATA  WILL  BE  OVERWRITTEN " 
INPUT  "ENTER  THE  NAME  OF  THE  DATA  FILE  TO  CREATE",  FILES 
IF  FILES  = ""  THEN  GOTO  FIT  FIT 
OPENIT: 

OPEN  FILES  FOR  OUTPUT  AS  #3 
PRINT  #3,  DATES 
PRINT  #3,  TIMES 
PRINT  #3,  "Pump  Set  ="; 

PRINT  #3,  SET 
PRINT  #3,  "Elec.Power  ="; 

PRINT  #3,  USING  "###.###";  P 
PRINT  #3,  "Act.Power  ="; 

PRINT  #3,  USING  "###.###";  PTOT 
PRINT  #3,  "Lost  ="; 

PRINT  #3,  USING  "###.###";  DELP 
PRINT  #3,  "Flowrate  =  "; 

PRINT  #3,  USING  "###.####";  m 
PRINT  #3,  "Density  =  "; 

PRINT  #3,  USING  "###.####";  den 
PRINT  #3,  "Viscosity  = "; 

PRINT  #3,  USING  "###.####";  vis 
PRINT  #3,  "Conductivity  = "; 

PRINT  #3,  USING  "###.####";  kf 
PRlNT#3,"Prandtl  =  "; 

PRINT  #3,  USING  "###.####";  Pr 
PRINT  #3,  "Cp="; 
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PRINT  #3,  USING  "###.####";  Cp 
PRINT  #3,  "Velocity  =”; 

PRINT  #3,  USING  "###.####";  VEL 
PRINT  #3,  "Reynolds  ="; 

PRINT  #3,  USING  "###.####";  Re 
PRINT  #3,  "Flux  ="; 

PRINT  #3,  USING  "####.###";  FLUX 
PRINT  #3, "  Nusselt  ="; 

PRINT  m,  USING  "###.###";  Nu 
PRINT  #3,  "Colburn  j="; 

PRINT  #3,  USING  "###.####";  Je 
PRINT#3,"Vflmtr="; 

PRINT#3,V(0) 

PRINT#3,"Vheat="; 

PRINT#3,V(1) 

PRINT  #3,  "Vresis="; 

PRINT#3.V(2) 

PRINT  #3,  "’llennocouples(C)" 

FORI  =  OTO  17 

PRINT  #3,  (I  +  100);  TMP(I);  (I  +  400);  TEMP(I) 

NEXTI 

PRINT#3,"Tsurf="; 

PRINT  #3,  USING  "###.##";  TSURF 
PRlNT#3,"Tfilm="; 

PRINT  #3,  USING  "##.###";  FILMTEMP 
PRINT  #3,  "Tin  ="; 

PRINT  #3,  USING  "##.###";  INLET 
PRINT  #3,  "Tout  ="; 

PRINT  #3,  USING  "##.###";  OUTLET 
PRINT  #3,  "Delta  =  "; 

PRINT  #3,  USING  "##.###";  DELTA 
PRINT#3,"Tdiff= "; 

PRINT  #3,  USING  "##.###";  TEMDIF 
PRINT#3,"Tcover="; 

PRINT  #3,  USING  "##.###";  TCOVER 
PRINT  #3,  "Tinsulation  = "; 

PRINT  #3,  USING  "##.###";  TINSULATION 
QUES: 

INPUT  "DO  YOU  WISH  TO  TRY  ANOTHER  RUN  (Y/N)  ?  ",AGAIN$ 

IF  AGAIN$  =  "  Y"  THEN  GOTO  STARTIT 

IF  AGAINS  =  "y"  THEN  GOTO  STARTIT 

IF  AGAINS  =  "N"  THEN  END 

IF  AGAINS  =  "n"  THEN  END 

GOTO  QUES 

END 
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APPENDIX  C.  PROPERTIES  OF  BRAYCO  MICRONIC  889  (PAO) 


The  properties  of  Brayco  Micronic  889  Polyalphaolefin  taken  from  Ref  10  are  tabulated 
below.  The  thermal  conductivity  of  Brayco  Micronic  889  (PAO)  was  0.083-0.085  BTUIhr.ft'^ 
during  the  heat  transfer  experiments. 


TEST  (ASTM) 

DESCRIPTION 

RESULT 

D-445 

1  Kineaiatic 

Viscosity 

CID^S  m-^2/s) 

@  100 'C 

1.7 

0  40  ®C 

52 

0  -40<’C 

260 

0  -54  “C 

1200 

D^287 

Density 

(g/jal) 

0  O^C 

0.811 

0  20 'c 

0.794 

0  40  "C 

0.777 

0  100 

0.723 

0  160 “C 

0.  661 

D-2766 

Specific  Heat 

(cal/g  C) 

0  -18»C 

0.  49 

0  10 

0.52 

0  38“C 

0.54 

0  93'C 

0.58 

Table  9.  Properties  of  Brayco  Micronic  889 [Ref.  10]. 
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APPENDIX  D.  SAMPLE  CALCULATIONS 


The  following  calculation  is  for  the  following  setpoint. 

Coolant :  Water. 

Electrical  Power  Input ,  P  Et- 100. 172  Watts. 

Flow  (%)  =  40. 

Mass  Flow  Rate,  m  =  0.0647  kg/s. 

Inlet  Temperature,  T1  =  20.05  °C 
Outlet  Temperature,  T2  =  20.26  °C 
Plate  Average  Temperature,  T3  =  24.02  '’C 

Temperature  Difference  Through  the  Insulation  Layer,  =  3.088  "C 
1.  Characteristic  Dimension 

Ac  =  sy.h 

Ac=  llJSmmx  15.24mm. 

^c=179.1mm2. 


Af=My.Ac 
4/=  19.5x179.1 
Af=  3492.5mm^. 

A=sxl+2xhxl+2xhxt 

^  =  (1 1.75  X  31.75)  +  (2  X  15.24  x  3 1.75)  +  (2  x  15.24  x  1.52) 
^=  1387.1mm2. 

n  -Ml  4x179.1  _ 

~  (A/I)  ~  (1387. 1/3 1 .75)  ^ 

2.  Properties  of  Water 

71+73 

Film  Temperature,  7/= — 2 — 

Using  Table  A.6  [Ref  9]; 

p=  1000.11  kg/m^ 
kf=  0.6005  W/m.K 

Pr=  6.625 

Cp  =  4180.51  J/kg.K 
p  =  9.623x10-^  N.s/m^ 


20.05+24.02 


=  22.04  “C 
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P  ”  1000.11 


=  6.469  X  10"^  mVs. 


_  Q  6.469x10-5  .  , 

,  V = ■T’  =  TTzrrr^  ~  ooiss  m/s. 
Af  3492.5x10-^ 


pxVxDf,  (1000.11)x(0.0185)x(16.39x10-5) 

=“i: - ^ 


Lost/Gamed  Heat  AP  = 


^insul^^  jfjsul 


where  A:*Ba/  (0.04  W/m.K),  is  the  thermal  conductivity  and  (0.01905),  is  the  thickness  of  the 


-(0.04)x(0.07742)x(3.088) 

^  ~  0.01905  ~  -0.502  Watts. 


P  =  Peic  —  AF 


F=  100.172  +  0.502  =  100.674  Watts. 
P  100.674 

Heat  Flux,  q”  =-j^  =  ^<^42  ~  W/m^ 


q”xD^  (1300.36)x(16.39xl0-3) 

“  kjxiI3-n}  ~  (0.6005)x(3.97)  “ 


Nuavg  8.94 

/■  =:—^ - =0.0151 

RenyPr^'^  f315.65)x(6.625)i^ 
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APPENDIX  E.  UNCERTAINTY  ANALYSIS 


The  accuracy  of  the  results  of  this  study  was  estimated  by  performing  an  uncertainty 
analysis.  The  theorem  established  by  S  .J.  Kline  and  F.  A.  McClintock  [Ref  11]  was  used  to 

calculate  the  uncertainties.  In  this  theorem  a  calculated  result  y  is  a  function  of  several 

independent  measured  variables,  {xi ,  xa, . ,x„ }.  Each  measured  value  has  some  uncertainty, 

{ui,  02, . ,Un},  and  these  uncertainties  lead  to  an  imcertainty  in  y,  which  is  called  %.  To 

estimate  %,  it's  assumed  that  each  uncertainty  is  small  enough  that  a  first-order  Taylor 
expansion  of  y(xi ,  xa , x„ )  provides  a  reasonable  approximation.  Under  this  approximation, 

y  is  a  linear  function  of  the  independent  variables  and  the  uncertainty  is  determined  by ; 


Here,  all  uncertainties  must  have  the  same  odds  and  must  be  independent  of  each  other 

/I  r  • 

For  y-A — j—  ,  the  uncertainty  takes  the  following  form; 

^3 

2  /  \2  /  \  211/2 


('^)  +(«^) 


The  uncertainty  analysis  performed  for  the  power  setting  of  100  Watts  and  the  Reynolds 
number  of  about  315.65  (for  water)  which  produced  a  Colburn  j  factor  of  0.0151  is  presented 
below: 

1.  Reynolds  Number  Uncertainty 


Re  = 


pxDfixV 

p 


where 


SRe 

Re 


I  f ; 
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For  Tfibn  =  22M°C 


p  =  1000. 1 1  kg/m^ ,  k  =  0.6005  W/m.K,  p  =  9.62  E-4  N-s/m^ ,  Pr  =  6.625 


Ac=0.0001791  m^. 
A/I  =  0.0437  m. 

Dh  =  0.01639  m. 

V  =  0.0 185  m/s. 


dAc  =0.00000962  m. 
8(A//)  =0.000707  m. 
6Z>/,  =0.000919  m. 
5F  =  0.0005  m/s 


The  imcertaimty  6Re/Re  =  0.0622  (6.22  %)  or  Re  =  315.65  ±19.63 


and 

j 

where, 

mu 

Nu 


For  ^”  =  1300.36  W/m^. 
vr(73-n)  =  3.97°c. 

=  8.94 
/  =  0.0151 

The  wncertainty  8j7j  =  0.093  or  9.3  % 


V’ =  43.32  W/m^. 

5vr=o.i  °c. 

=  0.625 
57  =  0.0014 


80 


LIST  OF  REFERENCES 


1 .  WING  AUNG,  "Cooling  Technology  for  Electronic  Equipment",  1988. 

2.  T.  Y.  TOM  LEE,  JAMES  A.  ANDREWS,  PETER  CHOW  and  DAVID  SAUMS, 
"Compact  Liquid  Cooling  System  for  Small  Movable  Electronic  Equipment",  TFKF. 
Transections  or  Components,  Hybrids  and  Manifacturing  Technology,  Vol.15,  No.5 
October  1992. 

3.  WEITING,  A.R.,  "Empirical  Correlations  for  Heat  Transfer  and  Flow  Friction 
Characteristics  of  Rectangular  Offset  Fin  Plate  Heat  Exchangers",  Trans.  ASME,  Journal 
of  Heat  Transfer,  Vol.  97,  Series  C,  488-490, 1975. 

4.  MANGLDv,  R.M.  and  BERGLES,  A.E,  "The  Thermal  Hydraulic  Design  of  the 
Rectangular  Offset  Strip  Fin  Compact  Heat  Exchanger",  Compact  Heat  Exchangers, 
Hemisphere,  123-140, 1940. 

5.  WEBB,  R.,  "Advances  in  Modeling  Enhanced  Heat  Transfer  Surfaces",  Heat 
Transfer  1994  Proceedings  of  the  Tenth  International  Heat  Transfer  Conference,  Vol.l, 
G.F.  HEWITT,  Ed.,  Taylor  Francis,  Bristol,  Pa.,  1994. 

6.  KAYS,  W.M.,  "Compact  Heat  Exchangers"  in  AGARD  Heat  Exchangers, 
AGARD-LS-57,  Advisory  Group  for  Aerospace  Research  and  Development,  NATO, 
Paris  1972. 

7.  S.V.  PATANKAR  and  C.  PRAKASH,  "An  Analysis  of  the  Effect  of  Plate 
Thickness  on  Laminar  Flow  and  Heat  Transfer  in  Interrupted-Plate  Passages",  Int.  J.  Heat 
Mass  Transfer,  Vol.24, 1801-1810, 1981. 

8.  LT.  JEFFREY  M.  MASTERSON,  "Heat  Transfer  Studies  on  a  Rectangular 
Channel  With  an  Offset  Plate  Fin  Array",  Thesis,  December  1993. 

9.  F.  P.  INCROPERA,  D.P.  DEWITT,  "Introduction  To  Heat  Transfer",  Second 
Edition,  Table  A-6.,  A22-A23, 1990. 

10.  CASTROL,  Speciality  Products  Division,  "Brayco  Micronic  889  Product  Data 
Sheet",  Irvine,  CA,  September  1994. 

11.  KLINE,  S.  J.,  and  F.  A.  McClintock,  Describing  Uncertainties  in  Single-Sample 
Experiments.  Mech.  Engr.  75;  3-8,  January  1953. 


81 


INITIAL  DISTRIBUnON  LIST 


1  Defense  Technical  Information  Center . 2 

Cameron  Station 

Alexandria,  Virginia  22304-6 1 45 

2  Library,  Code  52 .  2 

Naval  Postgraduate  School 

Monterey,  California  93943-5101 

3  Professor  M.D.  Kelleher .  3 

Department  Chairman,  Code  ME 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

4  Naval  Engineering  Curricular  Officer,  Code  34 .  1 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

5  Mr.  Tony  Buechler .  1 

Naval  Weapons  Support  Center 

Code  6042 
Crane,  Indiana  47522 

6  Ms.  Jennifer  Campbell .  1 

471  36*^  Ave. 

San  Francisco,  California  94121 

7  Deniz  Kuwetleri  Komutanligi . 4 

Personel  ve  Egitim  Daire  Baskanligi 

Bakanliklar  /  ANKARA  -  Turkey 

8  Istanbul  Teknik  Universitesi .  1 

Makina  Muhendisligi  Fakultesi 

Isi  Transferi  Kursusu 
Ayazaga  /  ISTANBUL  -  Turkey 


83 


Isi  Transferi  Kursusu 
ISTANBUL  -  Turkey 


10  Ortadogu  Teknik  Universitesi  .. 
Makina  Muhendisligi  Fakultesi 
Isi  Transferi  Kursusu 
ANKARA  -  Turkey 

1 1  Mr.  Fevzi  Komral  . 

Gencosman  Mah.  Guzeller  Sok. 
Ozen  Sitesi  B-BIok  No:  4 
Gungoren  ISTANBUL  -  Turkey 


